Cavity Light Bullets: 3D Localized Structures in a Nonlinear 

Optical Resonator 

Massimo Brambilla*, Lorenzo Columbo, Tommaso Maggipinto 
INFM, Dipartimento di Fisica Interateneo, 
Universita e Politecnico di Bari, via Orahona 4, 70126 Bari, Italy 

(Dated: today) 

Abstract 

We consider the paraxial model for a nonlinear resonator with a saturable absorber beyond the 
mean-field limit and develop a method to study the modulational instabilities leading to pattern 
formation in all three spatial dimensions. For achievable parametric domains we observe total 
radiation confinement and the formation of 3D localised bright structures. At difference from freely 
propagating light bullets, here the self-organization proceeds from the resonator feedback, combined 
with diffraction and nonlinearity. Such "cavity" light bullets can be independently excited and 
erased by appropriate pulses, and once created, they endlessly travel the cavity roundtrip. Also, 
the pulses can shift in the transverse direction, following external field gradients. 
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The competition between transverse diffraction and nonlinearities in optical systems (res- 
onators, systems with feedback mirror or counterpropagating beams) leads to the formation 
of cavity solitons (CS) appearing as bright pulses on a homogeneous background [l[Q]- The 
possibility of exploiting CS as self-organized micropixels in semiconductor vertical microres- 
onators has been recently proved after previous observations and predictions in various 
classes of macroscopic optical systems (see e.g. 0, l^l). Such phenomena occur when the 
mean-field-limit (MFL) is valid, so that in the propagation direction z (coinciding with the 
resonator axis) there is no spatial modulation of the field envelope. In this work we show 
how, beyond the MFL, an analogous formation of soliton-like structures confined in the 
transverse plane and in the longitudinal direction can be predicted. We call these structures 
Cavity Light Bullets (CLBs). They appear as bright stable pulses, spontaneously stemming 
and self-organizing from the modulational instability (MI) of a homogeneous field profile, 
and travelling along the resonator with a definite period. 

In the past, the temporal dynamics of a laser field, linked to the competition of longitu- 
dinal modes and without the MFL, has been extensively studied 0, Q. Those approaches 
included a plane wave approximation and could not account for pattern formation. On the 
other hand there are relatively few works dealing with pattern formation beyond the MFL 
in an optical resonator. An exact treatment has been proposed in |^ but negligible field 
absorption is considered, i.e. the field intensity is assumed constant along the 2;-axis. Other 
approaches proceed from a model proposed in P], where 2nd-order dispersion is the main 
mechanism for longitudinal self-modulation of propagating pulses , and the adaptation to 
resonator systems is heuristically provided with a formalism similar to MFL models . There, 
the formation of light bullets is shown jiol. Other theoretical studies have considered 
parametric resonators (x^ nonlinearities) where, though, the mechanism of structure confine- 



12| rather than pattern localization, and experimental 



ment is due to domain wall locking 
confirmations are still missing. 

On the other side, the self-collapse of a pulse, leading to formation of light bullets freely 
propagating in a medium, has been extensively studied by various authors [l^ . Moreover, in 
this case, the phenomenon does not involve a MI in the sense of Ref. lj| or the localization 
of patterns emerging thereof: the formation of one or more bullets involves a single pulse 
undergoing diffraction and dispersion. The CLBs presented in our work, conversely, are 
stable solutions of the intracavity field profile, superimposed to a nonmodulated background. 
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In this Letter, we analyze a well-established model for a resonator with a saturable 
absorber Q to which diffraction in the transverse plane is added U- The omission of 




the atomic dynamics, at difference from allows us to simplify the spatiotemporal field 
behavior avoiding competition between different timescales (field and atoms) 
After casting the model, we extend an approach developed in the 80s 
diffraction and analyze the MI domains and the pattern formation stemming thereof. We 
apply checkpoints using existing literature both close and far from MFL conditions. We 
then find conditions for radiation confinement and report about Cavity Light Bullets in a 
parametric domain not far (apart from the exclusion of the MFL) from [ij]. Finally, we 
show how CLBs can be excited and their motion controlled via field gradients. 

In the paraxial and Slowly Varying Envelope Approximation (SVEA), but without any 
hypothesis on the longitudinal profile of the intracavity field, the equations governing the 
spatiotemporal dynamics of a two level saturable absorber in a nonlinear unidirectional 
ring resonator driven by a coherent field, can be cast as done in and after adiabatic 
elimination of the atomic variables (good cavity limit) we have 

IdF OF -aF(l-iA) i „o ^ 

c dt dz l + A2 + |F|2^2fco ^ ^' 
(and complex conjugate) with the boundary condition 

F(x, y, 0, t) = TY,^, + RF{x, y, L, t)e-'''^ (2) 

where Ymj, F denote respectively the normalized envelope of the input beam (CW and 
transversely homogeneous) and of the intracavity field; A is the scaled atomic detuning, 60 
is the scaled cavity detuning, a is the absorption coefficient per unit length at resonance. 
T and R are the transmission and reflection coefficients, respectively, of the entrance and 
the exit mirror, being T + R = 1; x and y are the transverse Cartesian coordinates and 
z the longitudinal coordinate, being z = and z = L the input and the output mirrors, 
respectively. Here, for simplicity, we consider the equivalence between the cavity {C) and 
the medium's (L) lengths, i.e. the active medium completely fills the cavity. 

The transverse laplacian V^^ describes diffraction in transverse plane {x,y), while OF /dz 
accounts for the variation of the field envelope along the cavity axis. 

The stationary, transversely homogeneous solution is obtained by setting V^-F = and 
dF/dt = in Eq. ((TJ) and by solving the resulting ODE; it turns out that it can be 



multivalued. The stationary intracavity intensity |-Fst(-2)P is monotonically decreasing from 
input to output mirror 0, |2l|. Then, extending the approach of Lugiato et al. to include 
diffraction, we can study the stability of the solutions Fst{z) = pst{z)e''^''*^^\ {Pst,0st £ T^) 
against perturbations which are modulated along z and in the (x, y) plane; i.e. we consider 
a perturbed field F{x, y, z, t) as follows 

Fstiz) + 5Fix, y, z, t) = Fst{z) + 5/(z)e"*e*('=^^+'=-^) (3) 

with Sf{z) = e*^=* iSp{z) + ipstS9{z)), (to first order, with Sp,S9 G 71). and ky represent 
the components of the transverse modulation in the Fourier space. The different ansatz on 
the spatial dependence of the perturbations reflects the symmetry invariance in the [x, y) 
plane and the lack thereof along the axis. At difference from the MFL treatment, due the z 
dependence of pst and Ogt , we cannot obtain a simple algebraic equation for the eigenvalues 
A after linearizing Eq. (Q) around the imposed solution Fst{z) for the small perturbation 
6F{x,y, z,t). As it turns out, we must first find two linearly independent solutions ri{z), 
r2{z) and Ui{z), U2{z) for the following system of two coupled 1st order ODEs, as obtained 
from insertion of ansatz Q into Eq. (Q): 

dr [I l\ [ klA \ 

du /A klA \ f 1 \ 
dX='[A + Ak^)+''[2x) 



where X{z) = plt{z), r{z) = 5p{z)e'^, u{z) = 59{z)pstiz)e'^ , A = 1 + + X{z) and 
kj_ = k^ + ky. As the field must fulfill the constraints imposed by the boundary conditions, 
once such solutions are found, the perturbed field F{x,y, z,t) can be evaluated, and then 
substituted into Eq. (j2)). We rehandle the latter equation along the lines of and finally, 
we obtain a quadratic equation for 77 = i?e~^ 

7]^Wi + r]W2 + W3 = 0, (6) 

where the VTj's depend on 60, 9st{z) and on ri^2{z = 0), 7-12 (-2 = L), Mi 2(^2 = 0), Mi 2(^2 = L). 
The eigenvalues can thus be obtained as 

A± = (-^) (- ln(r/±) + ln(i?) - i2nn) , n e (7) 

where it is evident that the condition Re{X±) > 0, allows for a numerable set of unstable 
longitudinal modes to destabilize the stationary solution, for each {k^, ky) wavevector. From 



a physical view point we can then expect a complex unstable dynamical behavior of the 
system associated to a nonlinear competition among transverse and longitudinal modes. A 
detailed treatment of the stability analysis and of the guidelines for accessing a paramet- 
ric domain where localization in 3D could be achieved, fall beyond the letter format: an 
extended paper is currently in preparation. 
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We tested this stability analysis by reproducing the instability scenario predicted in 
when the values of A, aL and 5q approach the mean field limit (according to definition 
introduced in the number of longitudinal modes, setting off the instability of the 

stationary solutions, reduce to one. 

Beyond the MFL, we could use existing literature as a checkpoint, to reproduce some 
patterns previously reported far from the MFL, e.g. the multiconical character of the MI 
with a parametric set close to that considered in Patrascu et al. Q . 

The numerical integration of the dynamical equation Eq. has shown in this case the for- 
mation of three dimensional patterns strongly irregular in both space and time; this behavior 
agrees with the preliminary analyses suggesting chaotic dynamics and spatial decorrelation 
of patterns as the input field intensity increases beyond the MI threshold. 

In order to meet stable global patterns and 3D localization of light structures, the path 
we induced from our analyses was to limit the nonlinear radiation-matter coupling, without 
reducing absorption, and also remaining fairly far from the MFL. A condition we obtained 
by choosing: A = 0, aL = 1.2, T = 0.1, 6o = 0.1. In Fig. Q we report for this set of 
parameters the stationary transversely homogeneous state curve by plotting the intensity of 
intracavity field on the exit window / = versus the input filed amplitude, and the 

instability domain {Re{X±) > 0) obtained using the stability analysis just described. 

In this case, the simulations show, for example, a transverse intensity profile of isolated 
peaks, irregularly distributed in space and oscillating in time . Seen in 3D, we identify a 
series of structures, emerging from a subcritical bifurcation of the first positive branch of 
the steady-state curve (see Fig. Q). They are confined in the transverse plane but also have 
a limited, variable length in the z direction and travel the resonator with a period ~ ^ 
(see Fig. Efa)). This is thus a valid example of spontaneous self-organization of a complex 
optical system in 3D (and time). Moreover, we can observe a spontaneous reduction of the 
length of such structures, until a stable (and minimal) limit length is achieved (see FigI2fb)). 
These structures are what we call CLBs. 
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For input field values where structures coexists with the lower stable homogeneous branch, 
delimited in FigQ by the two arrows, by means of the usual technique adopted to switch 
on/off 2D Cavity Solitons in the mean field regime 22|, we added to the input field a suitably 
narrow and short Gaussian pulse, to locally realize a portion of the spatially modulated 
solution. 

By varying the Gaussian pulse duration, intensity and phase, we managed to "write" and 
"erase" the CLBs, i.e. structures confined in all three spatial dimensions (analogous to the 
ones reported in Fig. ^ which make a round-trip in a period ~ when the addressing 
pulse is much longer, the length of the structure reaches the full resonator's length and the 
stable structure emerging thereof becomes the 3D analog of the 2D cavity soliton. 

The computational most demanding simulations have been carried on in a simplified 
system with only one transverse dimension, x, after checking the coherence of the pattern 
scenario between the 3D and {z, x) systems. We show that is possible to create several inde- 
pendent CLBs using input Gaussian pulses aimed at any transverse position, provided they 
are separated at least by a critical distance on the order of the CLB diameter. Moreover, by 
superimposing to the holding beam a transverse phase modulation, we observe a slow lateral 
drift of CLBs climbing the intensity hill, induced by the input phase gradients themselves, 
towards their maxima (see Fig. 

We believe that the lower limit in the self-organised CLB length is linked to some char- 
acteristic length of the system, but we are still missing an unquestionable analytical result. 
We report that when longer pulses are applied in the input field, one seems to achieve longer 
CLBs (in the z direction); an organized report about the latter issues will be published 
elsewhere. The CLB reported in this work can be highly appealing for optical information 
processing, making possible the architecture of self-clocked, self-organised, reconfigurable 
pixels, which can encode all-optical information both serially (in the longitudinal trains of 
CLBs) and parallely (in the transversal arrays of CLBs). With respect to manipulation 
of 2D arrays of CS (proposed e.g. in [3, Q) one can figure out additional controls (e.g. 
transversely injected fields) to change the phase of CLB trains and thus manipulate infor- 
mation contents, or to increase the "input" channels, which can also be seen as logical gates' 
operands. 

On different grounds, CLBs lend themselves to similar applications as usual light bullets 
do, namely stroboscopes for atomic/molecular dynamics. Presently (because of the SVEA 
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in our model) our timescales are still too long for this, but the benefit unique to cavity- 
sustained structures is that they can be seen as a coherent state of the radiation in the 
cavity that may repeatedly interact with (e.g.) a quantum system (atom/condensate). In 
this case one could speak of a "quantum stroboscope" provided the coherences in the system 
are long compared to the cavity round-trip. 

This work was supported by the MIUR National Project "Formazione e controllo di 
solitoni di cavita in microrisonatori a semiconduttore" 
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FIG. 1: Right: Instability domain. Left: Steady state curve of the homogeneous solution and 
results of numerical simulations; continous and dashed lines refer to stable and unstable portion of 
the homogeneous branch, respectively. The ordinate of the asterisk and the solid triangles symbols 
corrisponds respectively to the mean values of the maximum intensity / of cavity light bullets and 
longitudinal rolls and filaments. /_ and /+ correspond to the minimum and the maximum intensity 
of the instability region and hence represent the ordinates of the limiting points A and B ond the 
stationary curve. The arrows delimit the interval of Yinj where the modulated solutions coexists 
with the stationary homogeneous one. 
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FIG. 2: Example of spontaneous self-organization of the system. Gray scale isosurface of intracavity 
field for Yinj = 6.74. 3-D patterns in the intensity profile (a) can collapse in localized structures (b) 
which make a round cavity trip in a period of ~ L/c. Note that in all numerical simulations the 
transverse spatial unit is given by the quantity L/2koT, while the longitudinal variable is scaled 
to L. Finally the time is scaled to the quantity L/cT so that the interval between the frames 
represented in this figure is around 545 time units. 
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FIG. 3: One transverse dimension case. Gray scale plots of intensity intracavity field profile; white 
represents the maximum values, (a, b, c) Sequence of frames showing the writing process of two 
indipendent CLBs (the whole sequence is 23 t.u. long); (d, e, f) CLB motion induced via a phase 
gradient in the holding beam (the whole sequence is 520 t.u. long). 
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